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摘要 ”利用 CHAMP 卫星 矢量 和 标量 地 磁 测 量 数据 ,通过 三 维 Taylor 多 项 式 模型 建立 了 2010. 0 年 中 国 及 邻近 地 
区 在 300 km 高 度 附近 的 X.Y.Z 和 下 分 量 的 三 维 磁场 模型 . 为 了 比较 验证 所 建 模型 ,分 析 了 其 与 Taylor 多 项 式 模 
型 (二 维 模型 ) 不 同 截断 阶 数 所 对 应 的 均 方 偏差 (RMSE) 、 残 差 及 分 布 等 .结果 表明 ,三 维 模型 每 一 阶 的 RMSE 和 残 
差 绝 对 值 的 平均 值 均 要 比 二 维 模型 的 小 约 415%. 由 于 采用 了 系数 完全 展开 的 建 模 方式 ,三 维 模 型 的 系数 数量 约 为 
二 维 模型 的 2 倍 ,三 维 模型 较 低 的 截断 阶 数 可 以 反映 更 多 的 地 磁 信 息 . 本 研究 中 的 5 阶 三 维 Taylor 模型 基本 可 达 
到 8 阶 Taylor 模型 的 精度 . 两 种 模型 绘制 的 地 磁场 及 残 差 分 布 有 较 好 的 一 致 性 . 
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Abstract The inherent shortcoming of most regional models is the lack of altitude information. 
In order to overcome this, the objective of this study is to establish the three-dimensional (3D) 
model of element X, Y, Z, and F over Chinese and its adjacent regions in 2010. 0 at the height of 
about 300 km. Based on the vector and scalar geomagnetic measuring data of CHAMP satellite, 
associated with 3D Taylor polynomial method. we try to create the 3D model over China and 
calculate the spatial distribution of geomagnetic field region. 3D model not only has the 
advantages like simple calculation and convenient to use, but also takes into account the height. 


In addition, the selection of truncation level of 3D model was also discussed. In order to verify 


基金 项 目 ”国家 自然 科学 基金 (41404053,41174165) ,公益 性 行业 (气象 ) 科 研 专项 (GYHY201306073) ,江苏 省 高 校 自然 科学 研究 基金 
(14KJB170012) 以 及 大 学 生 创 新 创业 训练 计划 项 目 基金 联合 资助 . 

作者 简介 蒋 勇 ,主要 从 事 数 值 模拟 及 空间 天 气 建 模 方面 的 研究 . E-mail:jianguk@163. com 

* 通讯 作者 ” 冯 疹 , 男 ,1980 年 生 , 博 士 ,讲师 ,人 研究 方向 :地 磁场 建 模 以 及 地 磁场 与 空间 天 气 的 关系 . E-mail: frank_feng8848@163. com 


3122 地 球 物 理学 报 (Chinese J. Geophys. ) 58 卷 


the constructed 3D model, a comparison between 3D and Taylor polynomial (2D) model was 


taken, we analyzed the Root-mean-square error (RMSE) and the residual distributions between 


3D and 2D models in different truncation level. 


Results showed RMSE and the mean absolute residuals of each level of the 3D model are 
about 45% smaller than that of 2D model. The 3D model of degree 5 basically has the same 


precision as that of 2D model of degree 8 in this study. Geomagnetic distributions and residuals of 


3D and 2D models have high consistency. After comparing with 2D model. the 3D model 


considers the altitude information, so the precision is really improved; the coefficients show that 


3D model of low degree can reflect more information and has the similar precision when compared 


to 2D model of high degree. 
Keywords 


1 Sle 


区 域 地 磁场 模型 用 以 反映 某 个 区 域 的 地 磁场 的 
时 空 分 布 特征 . 相 比 于 全 球 地 磁场 模型 (Hurwitzd 
et al., 1966; IAGA, 2003, 2005, 2010; 8 È &, 
201 4a) ,结合 高 密度 的 地 面 实测 数据 ,区 域 模型 能 
精确 地 反映 区 域内 的 地 磁场 分 布 , 如 地 磁场 在 特定 
区 域内 所 呈现 的 一 些 变化 规律 和 异常 分 布 . 

常用 的 区 域 地 磁场 模型 包括 Taylor 多 项 式 模 
型 .曲面 Spline 模型 ( 安 振 昌 和 徐 元 芳 ,1981; 安 振 
昌 等 ,1982)、 矩 谐 模 型 (Alldredge,1981)、Legendre 
模型 .自然 正 交 分 量 CNOC) 模 型 、 球 冠 谐 模 型 
(Haines,1985) 等 . 国内 外 学 者 已 利用 这 些 区 域 模型 


Geomagnetic model; 3D Taylor polynomial; CHAMP satellite; Root-mean-square error 


(1992) 运 用 Taylor 多 项 式 模 型 计算 了 1950—1985 
年 中 国 地 区 地 磁场 长 期 变化 模型 . 高 金田 等 (2005) 
利用 2003 年 中 国 地 区 的 地 磁 台 站 数据 建立 Taylor 
多 项 式 模型 对 中 国 地 区 的 地 磁场 进行 了 人 研究 . YE 
等 (2010c) 利 用 2000 年 中 国 地 磁 台 站 和 复 测 点 数 
据 , 根 据 Taylor 多 项 式 模型 计算 了 2000 年 中 国 地 
区 地 磁场 . 

尽管 存在 上 述 优点 ,但 Taylor 模型 是 二 维 模 
型 , 建 模 时 只 考虑 了 经 度 与 纬度 ,而 忽略 了 高 度 因 
K. 地 磁场 的 强度 随 着 高 度 变 化 而 变化 ,一 般 以 
20 nT/km 的 速率 减少 ( 安 振 昌 等 ,1991b), 因 此 二 维 
模型 所 计算 的 地 磁场 强度 会 存 一 定 的 误差 . 针对 此 
问题 , 柳 士 俊 等 (2011) 提 出 了 三 维 Taylor 多 项 式 模 
型 ,考虑 了 高 度 的 因素 ,为 区 域 模型 的 研究 提供 了 


做 了 很 多 研究 . 如 , Alldredge(1982) 利 用 和 矩 谐 模型 
人 研究 美国 东南 部 区 域 磁场 . Haines(1985) 利 用 球 冠 
谐 模 型 分 析 了 加 拿 大 区 域 1960 一 1983 年 的 磁场 变 
化 . Dizgit 等 (1997) 利 用 和 矩 谐 分 析 研 究 了 欧洲 区 域 
的 地 磁场 . Thébault 等 (2004,2006) (Thébault and 
Gaya-Piqu€é,2008) 将 原 有 的 球 冠 谐 模型 进行 了 改 
进 ,并 利用 改进 后 的 模型 研究 了 法 国 地 区 的 地 磁场 
(Thébault,2008). 国内 学 者 主要 针对 中 国 地 区 开展 
了 一 些 研究 . 徐 文 炮 和 朱 岗 昆 (1984) 利 用 和 矩 谐 分 析 
研究 了 我 国 及 邻近 地 区 的 地 磁场 . 高 金田 等 (2006 ) 
运用 曲面 Spline 模型 研究 了 1900 一 1936 年 中 国 区 
域 地 磁场 以 及 其 长 期 变化 的 分 布 . 陈 斌 等 (2011) 运 
用 球 冠 谐 模型 计算 2005 年 中 国 地 区 地 磁场 及 其 长 
期 变化 . 

相 比 于 其 他 的 区 域 模型 ,Taylor 多 项 式 模型 由 


种 新 的 方法 . 

在 以 往 的 研究 中 ,建立 和 计算 区 域 地 磁场 模型 
所 采用 的 数据 大 多 来 自 于 地 磁 复 测 点 、 台 站 数据 或 
者 利用 其 他 模型 所 模拟 的 数据 . 由 于 地 磁 复 测 点 和 
台 站 数量 有 限 , 且 空间 分 布 很 不 均匀 ,这 将 直接 影响 
所 建 模型 的 精度 . 另外 由 于 缺少 国界 外 测 点 ,使 得 在 
边界 区 域 上 的 模拟 会 出 现 异 常 , 即 边界 效应 . 一 般 解 
决 边界 效应 的 方法 有 添加 边界 外 区 域 的 补充 点 和 扩 
大 拟 合 区 范围 两 种 ( 杨 云 涛 等 ,2009). 补充 点 的 数据 
通常 来 自 于 其 他 模型 ,如 国际 地 磁 参 考场 (IGRF). 
这 种 数据 本 身 就 存在 一 定 的 误差 ,因此 会 进一步 影 
响 所 建立 模型 的 精度 . 扩大 拟 合 区 范围 后 ,虽然 在 拟 
合 区 中 心 的 高 精度 范围 扩大 ,但 拟 合 区 的 边界 异常 
仍 没 得 到 解决 . 

对 于 CHAMP 卫星 的 磁 测 数据 ,国内 外 学 者 做 


于 其 计算 简单 及 使 用 方便 ,得 到 了 广泛 的 应 用 . 如 安 
振 昌 等 (1991a) 利 用 Taylor 多 项 式 模 型 建立 了 
1950—1980 年 中 国 地 区 主 磁 场 模型 . 徐 元 芳 等 


了 大 量 的 研究 工作 . 如 Sabaka 等 (2004) 利 用 CHAMP 
等 卫星 数据 建立 了 第 四 代 地 磁场 模型 综合 模型 
CM4. Olsen 等 (2006) 利 用 CHAMP 等 卫星 数据 建 
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立 了 CHAOS 模型 系列 . Maus 等 (2007) 利 用 CHAMP 
卫星 数据 建立 了 岩石 圈 磁 场 模 型 (MF) 系 列 . BE E 
发 等 (2009, 2010) 基 于 CHAMP 卫星 数据 建立 的 
POMME-4. 2S 模型 研究 了 中 国 及 邻近 地 区 的 长 期 
变化 .长 期 加 速度 以 及 磁 异 常 的 分 布 特征 . 然而 ,将 
CHAMP 卫星 磁 测 数据 进行 中 国 大 陆地 区 的 地 磁 
场 研 究 至 今 还 鲜 有 报道 . 

针对 目前 研究 现状 ,认为 基于 CHAMP 卫星 所 
提供 的 高 精度 、 高 密度 的 全 球 地 磁 标 量 数据 和 矢量 
观测 数据 进行 中 国 及 邻近 地 区 的 地 磁场 三 维 建 模 对 
于 研究 中 国 , 乃 至 东亚 地 区 的 高 精度 地 磁场 形态 分 
布 具 有 重要 的 理论 和 应 用 价值 . 现 拟 首先 对 
CHAMP 卫星 数据 进行 数据 预 处 理 及 网 格 化 ,得 到 
中 国 地 区 卫星 地 磁 网 格 数据 ,然后 基于 不 同位 置 和 
高 度 的 卫星 数据 建立 中 国 及 邻近 地 区 的 三 维 
Taylor 多 项 式 模型 (三 维 模型 ) ,分 析 了 4 个 地 磁 要 
素 (X、Y、Z、F) 的 分 布 特点 ,研究 了 截断 阶 数 对 模型 
精度 的 影响 ,并 在 拟 合 值 . 残 差分 布 ,截断 阶 数 与 
RMSE 值 等 方面 与 常用 的 Taylor 多 项 式 模型 (二 维 
模型 ) 进 行 比较 分 析 . 


2 数据 与 方法 


2.1 数据 
地 磁场 为 一 向 量 场 , 地 球 空间 任意 一 点 都 包含 
7 个 地 磁场 要 素 ,分 别 为 地 磁场 总 强度 下 ,北向 分 量 
X, 东 向 分 量 了 ,垂直 分 量 Z, 水 平 强度 AL Rifa D 
All Tee FA I. CHAMP 卫星 提供 了 2001—2010 年 的 
地 磁 标 量 数据 (Ff) 和 矢量 数据 (X、Y、2Z), 所 有 数据 
由 德国 地 球 科学 研究 中 心 (GFZ) 网 站 (http:W 
www. gfz-potsdam. de/startseite/) 提供 . CHAMP 
卫星 于 2000 年 7 月 15 日 从 俄罗斯 的 Plesetst 卫星 
发 射 基 地 发 射 至 近 极 轨 、 环 形 、 高 约 454 km 的 轨 
道 . 卫星 搭载 的 Overhauser 磁力 计 (OVM) 测 量 标 
量 数据 ,人 磁 通 门 磁力 计 C(FGM) 测 量 矢 量 数据 ( 汉 郑 
等 ,2010a). 
由 于 所 测量 的 矢量 数据 是 在 FGM 参考 坐标 系 
下 获取 的 ,因此 还 需要 将 数据 转换 至 NEC( 北 向 、 东 
向 、 垂 直 ) 坐 标 系数 据 , 具 体 转 换 过 程 如 下 : 
Bygc E RvEce ITRE) i R arrer : Racrr-asc) 
e Reascercm 。 Brom > (1) 
其 Roascerom 为 将 FGM 传感器 参考 系数 据 转换 
至 卫星 光 具 座 (ASC) 参 考 系数 据 , 虽 然 角度 在 发 射 
前 已 经 确定 ,但 在 卫星 运行 过 程 中 会 不 断 变化 ,因此 


需要 进行 校正 ; Racrr-asc) 为 将 基于 level-2 的 卫星 
姿态 数据 由 ASC 参考 系数 据 转换 到 国际 天 球 参考 
系 (ICRF) 数 据 ; Ramrao 为 将 ICRF 参考 系数 据 
转换 到 国际 地 球 参 考 系 (ITRF) 数 据 ; Reece rer) 为 
将 ITRF 参考 系数 据 转 换 到 NEC 参考 系数 据 . 

为 了 尽量 减少 磁 扰 ,提高 模型 的 精度 ,本 文 参考 
常用 的 全 球 模 型 ,如 CM4 (Sabaka et al, 2004)， 
CHAOS3(Olsen et al,2010) 等 的 建 模 经 验 , 选 取 夜 
侧 卫 星 数据 , 即 18 :00 一 5 :00LT 期 间 数 据 . 并 从 中 
筛选 出 Kp 二 2,|D. | 二 20 nT 范围 内 的 数据 . 

本 文 拟 计算 并 分 析 2010.0 年 (2010 年 1 月 1 日 ) 
的 中 国 及 邻近 区 域 地 磁场 分 布 ,人 研究 范围 为 18 N— 
54°N,73°E—136°E, 为 了 获得 均匀 且 高 密度 分 布 的 
卫星 数据 ,通过 试验 ,发 现 对 于 某 一 时 间 点 前 后 60 
天 内 的 卫星 数据 经 归 算 后 可 获取 较为 理想 的 测 点 分 
布 ,然而 由 于 测 点 的 分 布 并 不 均匀 ,考虑 到 每 个 测 点 
在 三 维 建 模 中 的 权重 应 一 致 ,因此 对 归 算 后 的 所 有 
数据 进行 网 格 化 处 理 . 现 选取 2010. 0 年 作为 研究 时 
间 点 ,将 该 时 间 点 前 后 60 天 内 的 数据 进行 归 算 ,以 
经 纬度 1 为 间隔 ,在 每 个 格 点 处 取 其 上 下 左右 1 内 
所 有 卫星 数据 的 平均 值 进 行 网 格 化 . 对 剩余 较 少 的 
缺 测 格 点 , 则 进行 分 段 线性 插值 . 最 后 得 到 在 中 国 及 
邻近 地 区 的 卫星 网 格 值 , 共 为 2368 个 数据 点 . 2010 
年 为 CHAMP 卫星 观测 末期 ,卫星 轨道 高 度 较 初期 
有 所 下 降 ,2010.0 年 所 有 CHAMP 卫星 数据 的 高 度 
约 为 302. 59 一 312. 84 km. 

1 为 背景 网 格 数据 的 平面 分 布 图 ,每 个 数据 
格 点 都 均匀 分 布 .图 2 为 网 格 数据 的 三 维 分 布 ,从 低 
纬 到 高 纬 地 区 ,数据 点 高 度 呈 下 降 趋势 . 现 基于 这 些 
网 格 数据 建立 三 维 Taylor 模型 . 
2.2 Taylor 多 项 式 模型 

本 研究 基于 二 维 Taylor 多 项 式 模 型 建 Y 


< 
[1 
起 


T 了 T - 
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图 1 CHAMP 卫星 磁 测 数据 平面 网 格 图 
Fig.1 Two-dimenson distribution of CHAMP 


satellite magnetic grid data 
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图 2 CHAMP 卫星 磁场 数据 的 三 维 分 布 


Fig.2 Three-dimenson distribution of CHAMP satellite magnetic data 


Taylor 多 项 式 模型 . 二 维 模型 表达 式 如 下 : 


N i 
W= D Alepo Aaa, (2) 


i=1 j=1 


ERP OW 表示 任意 磁场 要 素 ; N APT pà 
为 各 地 磁 测 点 的 纬度 和 经 度 ; go ,ho 为 多 项 式 的 展 
开 原 点 , go = 36°N 5A, = 104. S°E; Ay 为 相应 各 分 量 
的 模型 系数 ,每 个 模型 有 N(N 十 1)/2 个 系数 ,所 有 
系数 通过 最 小 二 乘法 求 取 . 

(2) 式 的 系数 为 一 下 三 角 抢 阵 ,该 系数 展开 称 为 
经 典 展开 ,而 更 为 普遍 的 方法 为 进行 完整 展开 ( 方 阵 
展开 ) ,这 样 做 的 优点 是 可 以 较 低 的 截断 阶 数 反映 出 
更 多 的 地 磁场 信息 ,后 文 会 有 所 阐述 . 现 基于 二 维 模 
型 ,通过 添加 高 度 项 ,并 通过 系数 的 完整 展开 而 得 到 
三 维 Taylor 模型 ,其 表达 式 如 下 : 


W = x 5 S: (gp— p) T AS V hh), 
i=1 j=1 k=l 

(3) 
ERF W 表示 任意 磁场 要 素 ，N 为 截断 阶 数 ; p- 
hh 为 各 地 磁 测 点 的 纬度 、 经 度 和 高 度 ; po ,hu sho 为 多 
项 式 的 展开 原点 , m = 36 No = 104. 5°E, Ay = 
302.59 ~ 312. 84 km;Aj 为 模型 系数 ,每 个 模型 有 
N? 个 系数 ,所 有 系数 通过 最 小 二 乘法 求 取 . 


3 结 


3.1 截断 阶 数 的 选取 

Taylor 多 项 式 模型 为 基于 数据 点 的 拟 合 模型 ， 
截断 阶 数 的 选取 非常 重要 . 过 低 的 截断 阶 数 可 能 会 
导致 区 域 地 磁场 短波 长 信息 的 丢失 ,从 而 使 模型 无 


法 展示 高 精度 区 域 地 磁场 的 时 空 分 布 ,而 过 高 的 截 
断 阶 数 在 增加 计算 量 的 同时 ,还 会 使 结果 出 现 局 部 

通常 根据 如 下 两 条 依据 以 确定 合适 的 截断 阶 数 
( 安 振 昌 ,2001) : 

(1) 区 域 地 磁场 模型 的 均 方 偏差 要 小 于 IGRF 
的 均 方 偏差 ; 

(2) 随 着 截断 阶 数 的 增加 ,RMSE 逐渐 减 小 , 当 
RMSE 趋 于 稳定 时 ,可 确定 截断 阶 数 . 

RMSE 通过 下 式 计 算得 到 


5 CW moa es, W obs ) í 
RMSE = ， (4) 


7 一 上 
其 中 ,Wo 为 任意 地 人 磁 要 素 的 模型 值 , Wo 为 相同 
地 点 同一 地 磁 要 素 的 观测 值 , 为 观测 点 数量 . 
现 根 据 式 (4) 计 算 各 阶 三 维 模型 的 RMSE 值 ， 
结果 见 表 1. 
表 1 三 维 Taylor 多 项 式 模型 各 截断 阶 数 的 RMSE( 单 位 :nT) 
Table 1 The RMSE of each truncation order of 
3D Taylor polynomial model (units: nT) 


N X Y Z F 

1 5495 1453 10332 4710 
2 437.9 208 1108. 2 942.1 
3 39. 86 187. 88 101. 29 69. 27 
4 21.55 49,51 60. 16 45, 65 
5 15.72 15. 08 30. 51 17.53 
6 10445 492 14520 17917 
7 13657 773 18971 23353 
8 17225 945 23874 29406 
9 18946 1191 26371 32426 
10 21312 1285 29335 36263 
11 21824 1399 30127 37230 
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根据 表 1, 基于 CHAMP 卫星 数据 的 三 维 
Taylor 多 项 式 模型 的 RMSE 值 在 N<5 WY. X.Y, 
Z. FEM RMSE 随 N 的 增 大 而 迅速 减 小 ,平均 
下 降幅 度 达 到 了 65. 18%, 其 中 XX 分 量 由 1 阶 的 
5495 nT 减 小 为 5 WAY 15.72 nT,Y 分 量 从 1 阶 的 
1453 nT 减少 为 5 阶 的 15.08 nT, Z 分 量 从 10332 nT 
减少 为 30.51 nT, F 总 量 由 4710 nT 减 小 为 17. 53 nT; 
当 N>5 时 ,三 维 模型 的 RMSE 值 大 幅 增加 ,X 分 
量 为 10445 nT,Y 分 量 为 492 nT,Z 分量 为 14520 nT, F 
总 量 为 17917 nT. 这 是 由 于 阶 数 过 高 而 导致 的 区 域 
异常 现象 , 即 “ 龙 格 现 象 ”. 
根据 第 一 条 选取 依据 , 现 基 于 最 新 发 布 的 
IGRF12 计算 了 其 与 实测 点 数据 的 RMSE 值 ( 见 表 2). 


表 2 实测 点 与 ICRF12 的 RMSE 值 (单位 :nT) 
Table 2 The RMSE values between measuring points 
and IGRF12(units: nT) 


模型 X Y Z F 


IGRF12 4394. 12 595. 08 6175. 49 7595. 54 


根据 表 1 MA 2,4 22N<5 时 ,三 维 模型 的 
RMSE 值 小 于 IGRF12 的 RMSE 值 , 即 满足 选取 判 
据 (1). 综 上 ,三 维 模型 的 截断 阶 数 选 为 5. 


图 3 基于 5 阶 三 维 Taylor 多 项 式 模型 的 中 国 及 邻近 地 区 地 磁 分 布 (高 度 


3.2 中国 及 邻近 地 区 的 三 维 Taylor 多 项 式 模 型 的 

计算 和 分 析 

基于 不 同 高 度 的 CHAMP 卫星 网 格 点 ,利用 5 
阶 三 维 Taylor 多 项 式 模 型 计算 2010. 0 年 的 中 国 及 
邻近 地 区 地 磁场 网 格 值 . 模型 原点 为 研究 区 域 的 中 
心 , 即 po 二 36"N,N = 104, 5°E. 基于 所 得 系数 绘制 
的 卫星 高 度 处 XYZ F 四 个 要 素 的 三 维 分 布 见 图 3. 

根据 图 3,X 分 量 的 强度 随 着 纬度 的 升 高 而 逐 
渐 减 小 ,由 34140. 77 nT(18"N,59"E) 左 右 减 少 为 
14823. 98 nT(54°N,73°E) AA. ER E ly AKG 
纬度 方向 一 致 .Y Op et WY oD A WR BAR. KB oe 型 
分 布 ,强度 自 西 向 东 逐 渐 减 小 ,最 大 值 为 西北 部 (54 N， 
73°E) Mb AY 2449. 25 nT 左右 ,最 小 值 在 东北 部 (49"N， 
135°E) Mb) — 3527. 66 nT AA. Z 分 量 的 等 值 线 走 
向 基本 与 纬度 方向 一 致 ,强度 随 着 纬度 的 升 高 而 增 
大 ,由 低 纬 度 12349. 70 nT(18"N,136"E) 左 右上 升 
到 高 纬度 49657. 99 nT(54"N,99"E) 左 右 . 地 磁场 下 
总 量 的 分 布 与 2 分 量 类 似 , 强 度 随 着 纬度 的 升 高 而 
增 大 ,最 大 值 为 51884. 40 nT(54"N,104°E) 左 右 ,最 
小 值 为 33647. 11 nT(18"N,136"E) 左 右 . 由 于 采用 
的 是 卫星 网 格 数据 ,因此 不 存在 边界 效应 . 

为 了 比较 验证 所 建 模型 ,基于 相同 的 CHAMP 


100 11 
RIC) 


10°59 80 90 


1302. 59~312. 84 km) ;单位 :nT 


Fig. 3 Geomagnetic distribution of China and its adjacent regions based on 3D Taylor polynomial model 


with degree 5(altitude; 302.59~312. 84 km); units; nT 
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卫星 数据 ,对 所 有 高 度 都 取 平均 值 , 即 ho =307. 69 km. 
通过 输入 每 个 测 点 在 ho 处 的 经 纬度 及 X YY、Z 和 下 
总 量 的 强度 值 ,利用 二 维 Taylor 模型 计算 并 绘制 了 
中 国 及 邻近 地 区 的 地 磁场 分 布 (未 列 出 ). 当 N> 


58 4 


拟 合 值 要 略 好 于 二 维 模型 拟 合 值 . 通过 比较 ,可 以 发 
现 每 个 实测 点 的 二 维 以 及 三 维 模型 值 与 对 应 实测 数 
据 的 平均 误差 为 1%. 

然而 ,二 维 模型 是 一 种 平面 方法 ,所 有 数据 点 都 


时 ,二 维 模型 的 RMSE 值 逐 渐 趋 于 平稳 ,下 降 并 不 


要 求 在 同一 水 平 高 度 上 . 因此 二 维 模型 所 使 用 的 平 


明显 . 考虑 到 截断 阶 数 越 高 ,计算 量 越 大 , 且 两 种 模 
型 在 相同 阶 数 时 反映 的 信息 量 不 同 . 通过 RMSE 值 
的 比较 , 现 认 为 5 阶 的 三 维 Taylor 模型 与 8 阶 的 二 
维 Taylor 模型 都 能 达到 满意 的 模拟 精度 ,其 精度 接 
近 ( 后 文 有 详细 阐述 ). 比较 发 现 两 种 模型 所 得 到 的 
4 要素 地 磁场 无 论 在 分 布 还 是 强度 上 均 非 常 相似 ， 
故 不 再 将 二 维 模型 的 图 形 列 出 . 

为 了 比较 所 建 模 型 的 拟 合 效 果 , 现 将 部 分 测 点 
的 实测 值 及 其 8 阶 二 维 模 型 和 5 阶 三 维 模型 值 进行 
对 比 ( 见 图 4). 

根据 图 4, 两 种 模型 的 拟 合 值 与 实测 值 非常 接 
近 , 图 中 两 条 线 基 本 重合 ,而 Y 分 量 则 存在 一 些 差 
异 , 从 与 实测 值 的 接近 程度 及 变化 趋势 看 ,三 维 模 型 


4 
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均 高 度 (307. 69 km) 与 本 文 数据 (302. 59 一 312. 84 km) 
相 比 就 会 存在 士 5 km 左右 的 高 度 差 ,由 此 而 产生 的 
误差 是 不 可 避免 的 ,并 且 对 于 二 维 模型 是 无 法 计算 
这 部 分 误差 的 . 而 三 维 模型 的 结果 不 仅 考虑 了 高 度 
的 因素 而 且 能 与 实测 值 很 好 的 拟 合 . 因此 总 体 而 言 ， 
三 维 模型 的 拟 合 效果 是 令 人 满意 的 . 

为 了 进一步 地 观察 两 模型 的 差异 , 现 绘制 了 两 
种 模型 相对 于 实测 值 的 残 差 分 布 图 ( 见 图 5). 

通过 比较 发 现 ,两 种 模型 绘制 的 残 差 网 格 图 无 
论 在 分 布 以 及 强度 上 都 很 接近 . 并 在 各 分 量 的 一 些 
正 负 残 差 分 布 位 置 都 较为 相似 ,如 X 分 量 在 塔里木 
盆地 的 约 30 nT 的 正 残 差 分 布 和 内 蒙古 中 部 的 约 
一 40 nT 的 负 残 差分 布 .Y 分 量 在 中 国 地 区 东经 110” 
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图 4 三 维 模型 和 二 维 模型 的 模拟 值 与 实测 值 的 比较 ;单位 :nT 


Fig. 4 


The comparison among observed values and their modelling values based on 3D and 2D model; units: nT 
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图 5 三 维 模型 ( 左 ) 和 二 维 模型 ( 右 ) 的 残 差 分 布 图 ;(a 一 d) X,Y,Z, F;(e—h) X,Y,Z, 下 ;单位 :nT 
Fig.5 The residual distribution of 3D model (left) and 2D model (right); (a—d) X,Y,Z,F;(e—h) X,Y,Z,F; units; nT 


附近 约 一 33 nT 的 负 残 差分 布 .2 分 量 在 广西 南部 
约 一 70 nT 的 负 残 差 分 布 以 及 下 总 场 量 在 漠河 一 带 
约 40 nT 的 正 残 差 分 布 都 很 相似 . 

以 上 两 种 分 布 的 比较 很 难 发 现 明显 的 区 别 ,为 
了 更 为 直观 地 观察 两 种 模型 的 拟 合 效果 ,计算 了 两 
模型 相对 于 所 有 CHAMP 卫星 网 格 数据 的 残 差 . 比 
较 发 现 两 模型 的 残 差 在 变化 幅度 上 相差 不 大 ,但 对 
于 X 分 量 和 下 总 量 , 三 维 模型 要 明显 优 于 二 维 模 
型 . 三 维 模型 的 残 差 整体 要 小 于 二 维 模型 ,如 三 维 模 
型 X 分 量 的 平均 相对 误差 为 0.051%% ,要 小 于 二 维 
模型 的 0. 054%% ,三 维 模型 下 总 量 的 平均 相对 误差 为 


0.031% ,要 小 于 二 维 模型 的 0.043%, FRET 0.012%, 

为 了 进一步 说 明 两 种 模型 的 差异 ,通过 下 式 计 
算 两 种 模型 N 为 1 一 10 时 的 残 差 绝对 值 的 平均 值 
变化 ， 


W no — W os 
L Ža [Wawi aR (5) 
n 


ERE Ra 为 所 有 点 残 差 绝对 值 的 平均 值 , 其 他 各 
变量 的 意义 与 (4) 式 相同 . 

6 列 出 了 两 模型 在 N 为 1 一 5 时 的 残 差 绝对 
值 求 平均 的 结果 . 根据 图 6 可 知 , 随 着 N 的 增加 ,三 
维 模型 的 平均 残 差 的 下 降 速度 要 略 快 于 二 维 模 型 


Ra 
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平均 绝对 残 差 mT 


平均 绝对 残 差 mnT 


2 3 4 
截断 阶 数 N 
图 6 三 维 模型 和 二 维 模型 的 残 差 绝对 值 的 平均 值 变化 
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Fig.6 Mean values variation of absolute residual of elements X, Y, Z, F based on 3D and 2D model; units; nT 


的 ,在 N 相同 的 情况 下 ,三 维 模型 的 平均 残 差 要 比 
二 维 模型 平均 小 45. 09% ,这 说 明了 相同 阶 数 三 维 
模型 的 拟 合 精度 要 好 于 同 阶 的 二 维 模 型 

下 面 对 于 两 模型 的 RMSE 进行 分 析 . 首先 列 出 
了 二 维 模型 在 六 为 1 一 11 时 的 RMSE 变化 ( 见 表 3). 


表 3 二 维 Taylor 多 项 式 模型 各 截断 阶 数 的 RMSE( 单 位 :nT) 
Table 3 The RMSE of each truncation order of 2D Taylor 


— 


polynomial model (units: nT) 


N X Y Z PF 

1 5495 1453 10332 4710 
2 755 616.3 1677. 2 976.1 
3 309. 21 227. 27 302. 37 231. 63 
4 45. 28 69. 48 133. 87 82.7 
5 29.01 32. 49 38. 94 31. 37 
6 7.85 15. 43 32.1 25. 23 
7 7.23 13.51 29. 68 24. 86 
8 6. 96 12 29. 35 24.61 
9 6.77 11.81 29.19 24. 55 
10 6.71 11. 48 29. 06 24. 41 
11 4843 243 6814 8320 


根据 表 3, 基 于 CHAMP 卫星 数据 的 二 维 
Taylor 模型 在 截断 阶 数 N<6 时 ,XX、Y、Z\ 下 分 量 的 
RMSE 随 NN 的 增 大 而 迅速 减 小 ,平均 下 降幅 度 达 到 
了 60. 62% ,其 中 XX 分 量 由 1 BPAY 5495 nT WA 


6 阶 的 17.85 nT,Y 分 量 从 1 BAY 1453 nT 减少 为 6 
阶 的 15. 43 nT,Z 分 量 从 10332 nT 减少 为 32.10 nT, 下 


4 Sr 


Joh HE H 


H 4710 nT 减 小 为 25. 23 nT; 4 N 为 6 一 10 


时 ,RMSE 平均 减少 幅度 为 2. 96% ,逐渐 趋 于 平稳 ; 
当 N>10 时 ,二 维 模型 也 出 现 了 “ 龙 格 现象 ”, RMSE 
值 大 幅 增加 . 
为 了 进一步 比较 ,绘制 了 两 模型 当 N=1~11 
时 的 RMSE 变化 
5 时 ,三 维 模型 的 RMSE 值 均 小 于 二 维 模型 , 当 


N=5 时 ,三 维 模 


LA 7). 由 图 7 可知, 当 为 1 一 


型 中 X 分 量 和 下 总 量 的 RMSE fË 


都 要 远 小 于 8 阶 二 维 模型 的 X 分 量 和 下 总 量 , 而 2 


分 量 与 7 一 10 阶 


的 5 


有 相近 的 精度 ， 


的 二 维 模 型 值 非常 接近 ,三 维 模型 


阶 Y 分 量 则 略 大 于 8 阶 的 二 维 模 型 . 综合 考虑 
磁场 分 布 . 残 差分 布 以 及 RMSE 的 比较 , 现 认 为 5 
阶 的 三 维 Taylor 模型 与 8 阶 的 二 维 Taylor 模型 具 


精度 甚至 更 高 . 
通过 比较 发 现 ,三 维 模型 在 N 为 1 一 5 时 
RMSE 值 比 对 应 的 二 维 模型 平均 要 小 45.8%, 男 外 
根据 两 种 模型 的 系数 个 数 , 认 为 在 基于 CHAMP 卫 
星 数据 以 模拟 中 国 地 区 地 磁场 的 情况 下 , 较 低 阶 的 
模型 可 反映 更 多 的 地 磁场 信息 ,其 精度 相当 于 


三 维 


高 阶 的 二 维 模型 . 
为 了 进一步 观察 两 模型 ,并 考虑 到 存在 的 差异 ， 
列 出 了 N=3 时 两 模型 的 系数 , 见 表 4, 表 5. 


AYE X 分 量 和 下 总 量 上 三 维 模型 的 
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图 7 二 维 及 三 维 Taylor 多 项 式 模 型 的 RMSE 随 截断 阶 数 N 的 变化 ;单位 :nT 
Fig.7 The relationship between RMSE and N of elements X,Y,Z,F in 2D and 3D Taylor polynomial model; units; nT 
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表 4 3 阶 三 维 Taylor 多 项 式 模型 系数 
Table 4 The coefficients of 3D Taylor polynomial model while N is 3 


coeff X Y Z F coeff X Y Z F 

A 2. 6356 X 10! —1. 1184 10° 3. 6920 X 10! 4. 5378 X 10! A223 4. 0196X107? 9. 2792X107? —2. 2684X107? —1. 8737X107? 
Au 1. 8244 < 10° 2.836010! 一 3.0513 义 10! —4. 7863 10! A231 9. 7056X107? —1. 4293X107? —5. 9247X107? —1. 7627X107? 
Aus 9, 3890 X 10° 2.569910! —5. 4618X100 —5.7399X 10° Ag32 —2. 3759X1073 —3.7049X107? —2. 5003X107? —2. 7506X107? 
Aiz 8. 5038 X 10° —7. 708410! —5. 875610! —4. 0083 10! A233 2984X T03: -—1. 1612X 10 4. 2343 10-4 7. 4171X1074 
Aiz —2. 4694X107! 5. 7523X107! —1. 1463X100 —7. 0559X107! || Azı —5. 0897X 10° 2. 0020X 10° —1.4831X10! —3. 1636X 10° 
A2 —2.9190X10~! —1. 9139X 10° 7. 8564X107! 4. 7338X107! || As —1.185210° —6. 4037X1078 7. 7148X107! 1. 8777X 10° 
Ası —8. 4217X107! 1. 1318X107! —3. 4192 10° —3. 1678X 10° Aa3 —6. 4995X1074 —3. 3960X107? 2. 3336X107? 3. 8016X107? 
Ais2 9.0272X 1074 —1. 4263X107! 2.8775X10 . 3662X107? || Asn —1. 5394X107? 3. 6781X107? . 0011X107! 4, 6874X107? 
Als —4. 9521X1073 —8. 6342X107? —6. 7235X107? —7. 8745X107? || As . 5478X107? —2. 6794X1073 —1. 2601X1078 —2. 4206X107 
Asi —5. 6348X 10? —3. 2395X107! 9. 8169X 10? 4. 6313 X 10? A323 3. 4452 x 1074 3. 2986X1074 一 3.5591 义 10 上 —3. 8670104 
A212 4. 0267 X 10° 9.9163 10° —8.5600 107! 7.5899 X 107! || Assi 5. 2130X1074 —3. 3949X1078 . 2816X1078 —5. 9806X1074 
Asis —2.7306X10° 1. 1916X107! 1. 7439X 10° 4. 1502X 10° A332 5. 0862X1074 —2. 3039X1074 —5. 6094X1074 —1. 7193X1075 
A221 2. 0462 X 10° —2. 9027X 10° 1. 4973X107! 5. 4465X1071 || As —8. 4805X1076 —3. 1064X1075 —4, 6306 10~° 5. 6562X 107° 
A222 4. 7735X107? —6. 4648X107! 2, 2293 107! 9. 0850 X 107? 


#5 3 阶 二 维 Taylor 多 项 式 模型 系数 
Table 5 The coefficients of 2D Taylor polynomial model while N is 3 


coeff X Y Z F coeff X Y: Z F 
An 2. 6334 X 104 —1.1095X10° 3.676810! 4, 5286 X 104 Az —5.3071X10° 1.0536X 10° —1.3394X10! —2. 3098X 10° 
Ai 5. 0955 X 10? —3. 9682X107! 9.5088 xX 10? 4, 2726 X 10? Ax 1.8126 10° —2. 8480 x 10° 2. 4282X107! 6, 2468x107! 


Ais 6.7776 X 10° —7, 120210! —4, 890310! —3, 4906 10! Azs = —7.6290X 107! 1.0246 107! —3.0339X 10° 2. 9828 X 10° 
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由 表 4 和 表 5 可 知 ,两 种 模型 在 第 一 项 系数 上 
非常 接近 ,而且 在 数值 上 比 其 他 项 系数 大 2 到 3 个 
数量 级 ,因此 两 种 模型 在 拟 合 结果 以 及 分 布 都 非常 


(3) 本 研究 显示 了 较 低 阶 数 的 三 维 模型 ,其 
RMSE 与 较 高 阶 数 的 二 维 模型 接近 ,通过 比较 两 种 
模型 不 同 截断 阶 数 对 应 的 RMSE 值 ` 以 及 磁场 和 残 


接近 . 相 比 于 二 维 模型 ,三 维 模型 的 系数 个 数 约 为 二 
维 模型 的 4 倍 , 地 磁 模 型 的 系数 能 反映 地 磁场 强度 
等 信息 ,如 IGRE 模型 中 各 阶 系数 能 反映 偶 极 子 场 
和 各 阶 非 偶 极 子 场 的 强度 信息 , 当 N>13, AWER 
谐 模型 的 系数 即 反 映 了 不 同 波长 所 对 应 的 岩石 圈 磁 
场 的 强度 信息 . Taylor 多 项 式 模 型 的 各 阶 系数 从 地 
球 物理 的 角度 而 言 , 应 包含 了 区 域 地 磁场 不 同 尺度 


差分 布 ,认为 三 维 模型 每 一 阶 的 RMSE 和 残 差 绝对 
值 的 平均 值 均 要 比 二 维 模型 的 小 约 45%. 
通过 计算 还 发 现 二 维 模型 的 边界 畸变 很 大 ,而 
三 维 模型 的 边界 效应 相对 较 小 .原因 在 于 ,在 完整 展 
开 的 三 维 模型 中 ,我 们 取 的 截断 阶 数 为 N 不 是 很 
大 ， 龙 格 现象 "尚未 得 以 发 展 ( 柳 士 俊 等 ,2011)， 
二 维 Taylor 模型 只 考虑 了 经 度 与 纬度 ,忽略 了 


的 强度 信息 ,但 由 于 其 为 区 域 模型 ,并 不 具备 正 交 
性 ,因此 是 很 难 从 功率 谱 的 角度 定量 地 分 析 不 同 尺 


高 度 的 因素 . 因此 二 维 模 型 由 于 本 身 的 限制 ,所 计算 
的 地 磁场 强度 会 存在 一 定 的 误差 . 针对 上 述 缺点 , 结 


度 的 区 域 地 磁场 与 不 同 阶 系 数 的 关系 , 故 三 维 模型 
具备 更 多 系数 的 具体 含义 ,也 是 以 后 工作 中 需要 考 
虑 的 问题 . 


4 结论 与 讨论 


本 文 基于 CHAMP 卫星 的 网 格 数据 运用 三 维 
Taylor 多 项 式 模型 计算 并 绘制 了 2010. 0 年 中 国 及 
邻近 地 区 的 地 磁场 分 布 ,并 通过 对 比 ,分 析 了 两 种 模 
型 的 残 差分 布 、 截 断 阶 数 与 RMSE 等 .结果 显示 5 
阶 的 三 维 Taylor 模型 与 8 阶 的 二 维 Taylor 模型 具 
有 相近 的 精度 . 通过 比较 两 种 模型 的 拟 合 值 . 残 差分 
布 ,截断 阶 数 与 RMSE 值 等 ,显示 三 维 模型 都 要 优 
于 同 阶 的 二 维 模 型 ,并 认为 较 低 阶 的 三 维 模型 可 以 
得 到 与 高 阶 二 维 模型 相近 ,甚至 更 好 的 拟 合 效 果 , 并 
得 出 以 下 结论 : 

(1) 相 较 于 其 他 二 维 模 型 ,如 二 维 Taylor 多 项 
式 模型 .曲面 Spline 等 ,三 维 Taylor 模型 考虑 了 高 
度 因 素 , 在 中 国 及 邻近 地 区 ,海拔 高 度 为 一 4 一 8 km 
左右 ,而 地 磁场 一 般 随 高 度 呈 20 nT/km 左右 的 速 
率 衰减 ,因此 三 维 Taylor 模型 相 较 于 二 维 模型 有 了 
本 质 的 提高 .本 研究 使 用 了 高 精度 .高 密度 的 
CHAMP 卫星 的 标量 和 矢量 数据 ,基于 此 建立 的 中 
国 地 区 三 维 模型 , 相 较 于 IGRF11、 甚 至 球 冠 谐 模 
型 ,不 仅 计 算 方 便 , 而 且 精 度 相 当 甚 至 更 高 . 

(2) =4E Taylor 多 项 式 模 型 为 完整 展开 的 拟 合 
模型 ,其 系数 数量 为 相同 截断 阶 数 下 的 二 维 模型 的 


RN 信 , 约 为 2N 倍 ,因此 认为 三 维 模型 在 较 低 的 


截断 阶 数 下 可 以 反映 出 更 多 的 地 磁场 信息 . 然而 由 
于 系数 数量 随 截 断 阶 数 的 增加 呈现 指数 式 上 升 ， 
此 三 维 模型 也 较 易 出 现 “ 龙 格 现象 ”. 


果 显 示 三 维 模型 可 以 更 准确 地 描述 地 磁场 的 空间 分 
布 .之 所 以 采用 完整 展开 ,是 因为 可 在 较 少 的 截断 阶 
数 计算 出 更 多 的 系数 ,从 而 提高 效率 . 

由 于 缺乏 2010.0 年 的 地 面 实 测 数据 ,因此 研究 
只 展示 了 卫星 高 度 处 的 地 磁场 分 布 ,今后 藻 可 获取 
地 面 及 海洋 地 区 ,其 至 航 磁 数 据 , 可 尝试 从 不 同 高 度 
联合 建立 三 维 Taylor 多 项 式 模 型 ,该 模型 的 空间 适 
用 范围 为 地 面 至 卫星 高 度 , 因 此 适合 进行 空间 磁场 
的 应 用 研究 . 车 结合 T04(Tsyganenko and Sitnov， 
2005) .CM4 等 模型 ,可 进一步 模拟 区 域外 源 场 以 及 
岩石 圈 人 磁场 的 形态 分 布 ( 汉 彦 等 ,2010b;2014b). 综 
上 所 述 , 三 维 Taylor 模型 具有 计算 简单 .使 用 方便 、 
且 具 有 相当 的 模拟 精度 , 适 于 进行 区 域 三 维 地 磁场 
的 研究 . 

致谢 感谢 空间 天 气 学 国家 重点 实验 室 开 放 课 题 以 
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